A simple classical mapping of the spin-polarized quantum electron gas: distribution 

functions and local-field corrections. 



m 



X3 
O 

o 



> 

in 
o 

On 
O 
ON 



o 



X 
S3 



M.W.C. Dharma-wardanaEl and F. Perrot 

National Research Council, Ottawa, Canada. KIA 0R6 
Centre d'Etudes de Bruyeres le Chdtel, P. O. Box 12, 91689 Bruyeres le Chdtel, France 

(24 May 1999) 

We use the now well known spin unpolarized exchange-correlation energy Exc of the uniform electron 
gas as the basic "many-body" input to determine the temperature Tq of a classical Coulomb fluid 
having the same correlation energy as the quantum system. It is shown that the spin-polarized pair 
distribution functions (SPDFs) of the classical fluid at Tq, obtained using the hyper-netted chain 
(HNC) equation are in excellent agreement with those of the T = quantum fluid obtained by 
quantum Monte Carlo (QMC) simulations. These methods are computationally simple and easily 
applied to problems which are currently outside the scope of QMC. Results are presented for the 
SPDFs and the local-field corrections to the rsponse functions of the electron fluid at zero and flnite 
temperatures. 

PACS Numbers: 05.30. Fk, 71.10. -^x, 71.45. Cm 



The uniform interacting electron gas (UEG) is the 
"text-book" many-body problem of Ferini liquids in 
metals, plasmas or doped semiconductors^! It provides 
a model exchange-correlation potential for the densitw 
functional theory of inhomogeneous electronic systems.S 
The random-phase-approximation (RPA) to the proper- 
ties of the UEG provides an important "baseline" which 
could be reached by many techniques like quantum- 
linked cluster expansions, Green's functions, or linearized 
equations of motion. RPA is a reasonable approxima- 
tion when the "coupling parameter" F = (potential en- 
ergy)/ (kinetic energy) of the Coulomb fluid is smaller 
than unity. The coupling parameter for the UEG at 
T=0 becomes identical with the mean sphere radius rs 
per electron, i.e., for the 3-D case, F = = {S/A-KnY^^. 
Here n is the electron number-density per atomic unit of 
volume. 

Improving the RPA becomes nontrivial if the results of 
the theory are expected to satisfy sum rules and provide 
physically realistic pair-distribution functions (PDFs), 
i.e., g(r), of the quantum system. Diagrammatic meth- 
ods look for resummations that conserve the sum rules. 
Ward identities etc. This line of research is followed ki 
the work of Hubbard, Vosko, Langreth, Geldart etpal.B. 
Calculations have been given by MacDonald et al.D for 
the self-energy, while Green et al., Ichimaru et al. and 
Richardson et al—and many others have studied the re- 
sponse functions Jj 

The equation of motion (EQM) method, when pushed 
beyond the RPA requires some "ansatz" to decouple or 
"close" the equations. Q Singwi et al. (SSTL) used a phys- 
ically motivated classical analogy to close the EQM by 
introducing the electron-electron PDF determined self- 
consistently. There are also attempts to directly fit the 
response functions or the PDF's to a form cOjUStraincd to 
satisfy the sum rules and such requirements Jj In spite of 
a large effort, the calculation of the electron-gas g{r) by 



STLS and other methods at arbitrary spin-polarization 
(C) and T is difficult. In fact, STLS gives negative 
PDFs for sufhciently large in the metallic range and 
does not satisfy the compressibility sum rule. Even 
the attempts to directly fit the PDFs do not guarantee 
positive-definite a g(r) for important regimes of rs and 

tB 

Another approach beyond RPA is to use a correlated 
wavefunction tJj — FD where is a correlation factor 
and D is a Slater determinant. The best F is deter- 
mined variationally. This approach leads to the Jackson- 
Feenberg enocav— functional which can be examined in 
several ways.EHt3 The QMC techniques also use such a 
tp and leadj-fi.g., to the variational Monte Carlo (VMC) 
method.EjtJ Thus Ceperley and Alder promded E^cirs), 
now available in several parametrized forms .Mr 1 The cor- 
relation energy given by various QMC methods agree 
with one another to about the same extent as with, say, 
the results from the Fermi-Hyper-Netted-Chain(FHNC) 
method.tll. Thus, while it is easy to get good corrlation 
energies, it is quite difficult to get other propertis like the 
PDFS and local-field corrections (LEG) at arbitrary C, T 
and in the metallic range of r^. 

In this letter we present a computationally simple 
method for calculating the PDFs and other derived quan- 
tities (e.g, response functions) of the UEG, given the un- 
polarized Ec from whatever source. Since PDFs of clas- 
sical fluids are easily obtainable via the Hyper-Netted- 
Chain (HNC) procedure|la we model the electron gas as 
a classical fluid containing two species (up, down spins). 
We ask for the temperature Tq at which the excess free 
energy of a Coulomb fluid, obeying the classical HNC- 
integral equation, becomes equal to the correlation en- 
ergy Ec of the UEG at the same density and at T=0. The 
suffix q in Tq signifies that this temperature parametrizes 
the quantum many-body interactions in the UEG. Using 
the temperature Tq at each we obtain gij(r), spin- 



1 



polarized correlation energies Ec{rs, C), local- field correc- 
tions etc. While many interesting applications are pos- 
sible, here we treat the 3-D spin-polarized UEG at zero 
and finite T. 

The physical motivation for our approach comes from 
density functional theory (DFT) where the interacting 
electron gas is replaced by a noninteracting gas of Kohn- 
Sham (KS) particles whose wavefunction is a simple de- 
terminant. Here the philosophy is very different to the 
Feenberg appma|Ch which uses a correlated determinent 
for fermions.EJE3 In DFT the many-body potential is re- 
placed by a single-particle KS potential. Since the natu- 
ral energy parameter of the clasical ensemble is the tem- 
perature, we look for a temperature representation of the 
interactions. 

Consider a fluid of mean density n containing two spin 
species with concentrations Xi = rii/n. In the following 
we deal with the physical temperature T of the UEG, 
while the temperature of the classical fluid T^f is 1//3. 
Since the leading dependence of the energy on temper- 
ature is quadratic, we assume that T'c/= + 
This is clearly valid for T = and for very high T. This 
assumption will not be probed more deeply in this letter 
where the main effort is for T=0. Our objective is to de- 
termine Tq, the temperature equivalent of the quantum 
correlations in the UEG at T=0. 

The pair-distribution functions for a classical fluid at 
an inverse temperature (3 can be written as 

5ij (r) = exp[f3(l)i-i (^) - (^) " Czj (r) + Bi-j (r)] (1) 

Here (r) is the pair-interaction potential between 
species For two electrons this is just the Coulomb 
potential Vcouif) if there is no exchange interaction. If 
the spins are parallel, the Pauli principle prevents them 
from occupying the same spatial orbital. Following the 
earlier work, notably by Lado,t3 wc introduce a "Pauli 
potential", 'P(r). Thus: 



+^cou(r) 



(2) 



The Pauli potential 'P{r) will be discussed with the PDFs 
of the non-interacting UGE, i.e., gfjir). The function 
h{r) = g{r) -1; it is related to the structure factor S{k) 
by a Fourier transform. The c(r) is the "direct correlation 
function (DCR)" of the Ornstein-Zernike (OZ) equations. 



hij{r) = Cy (r) -I- T^sUs J dr'hi 



,(|r-r'|K,(r) (3) 



The Bij (r) term in Eq. |^ is the "bridge" term arising from 
certain cluster interactions. If this is neglected Eqs. ^ 
^ form a closed set providing the HNC approximation 
to the PDF of a classical fluid. Various studies have 
clarified the rolcjof i?(r) and its treatment via "reference" 
HNC equations.Ej B{r) is important when the coupling 
constant P exceeds, say, 20. The range of P relevant to 
this work (e.g., P ~ 4.5 even for rg = 10 ) is such that 
the HNC-approximation holds. The HNC-approximation 



suffers from a compressibility inconsistency (CI), i.e., the 
excess compressibility calculated from the small- fc limit 
of the short-ranged part of c(k) does not agree with that 
obtained from the excess free energy. This CI can be 
corrected by choosing Bij (r) suitably. 

Consider the non-interacting system at temperature T. 
The parallel-spin PDF, i.e, gii{r,T), will be denoted by 
g^{r) for simplicity, since g^j{r,T), i 7^ j is unity. Con- 
sider the paramagnetic case Xi = 0.5. Denoting (ri — Y2) 
by r, it is easy to show that 

g^Av) = |jSk„k.n(fci)r.(fc2)[l - e^C^-''^)-] (4) 

Here n{k) is the Fermi occupation number at the tem- 
perature T. Eq. ^ reduces to: 



g%{r) = l-Fl{r) 
Frir) = (6^V4) / n(fc) 



sin{kr) kdk 
■ 2^' 



(5) 
(6) 



The Fermi momentum is denoted by kp. Thus (?j'(r) 
can be obtained from the Fourier transform of the Fermi 
function. Then c"(r) can be evaluated from g^{r) us- 
ing the OZ|j:£lations. The T=0 case can be evaluated 
analytically^ 

Assuming that g^^ (r) can be modeled by an HNC fluid 
with the pair interaction l3V{r), dropping the indices, we 
have: 



g°{r) = exp[-l3V{r) + h°{r) - c°(r)] 



(7) 



The k-space direct correlation function c"(fc) at T=0 de- 
cays as Akp /Zk for small fc, showing that the r-space form 
c°(r) is long ranged. The "Pauli potential" 7'(r) is given 
by 



l3V{r) = -log[g\r)] + h^{r) - c\r) 



(8) 



We can determine only the product (3V{r). The classical 
fluid "temperature" 1//3 is still undefined and is not the 
thermodynamic temperature T . The Pauli potential is 
a universal function of rkp at T=0. Plots of (3V{r) and 
related functions are given in Fig. 1. It is long ranged 
and mimics the exclusion effects of Fermi statistics. As 
T is increased, the Pauli potential restricts to about a 
thermal wavelength and becomes more hard-sphere like. 

The next step is to use the full pair-potential ipij (r) , 
Eq. H, and solve the coupled HNC and OZ equations for 
the binary (up and down spins) interacting fluid. For the 
paramagnetic case, = n/2, we have: 



^ g-/9T'(r)5.,+yeou(r)+ft,j(r)-c,,(r) 
FT 



(9) 



h,M) ^ h^jir) (10) 
hii{q) = cii(q) + {n/2)[cii{q)hii{q) + Ci2{q)h2i{q)] (11) 
hi2{q) = ci2(q) + {n/2)[cii{q)hi2{q) + Ci2{q)h22{q)] (12) 

The Coulomb potential Vcouir) needs some discussion. 
For two point-charge electrons this is l/r. However, de- 
pending on the temperature T, an electron is localized 
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to within a thermal wavelength. Thus, following Morita, 
and Minoo et al.J^ we use a "diffraction corrected" form: 

Vcou{r) = (l/r)[l - e-'-/^*'-]; A*,, = {2^rnT,j). (13) 

Here m is the reduced mass of the electron pair, i,e, 
m*{rg)/2 a.u., where m*{rs), the effective electron mass. 
It is weakly dependent, e.g, ~0.96 for = 1. In this 
work we take r7i* = l. The "diffraction correction" ensures 
the correct behaviour of (712 (?" ~ 0) for all r^. 

In solving the above equations for a given and at 
r=0, we have Tcf=Tq. A trial Tq is adjust to obtain an 
Ec{Tq) equal to the known paramagnetic Ec{rs) at each 

Ts- 

Ec{Tq)^j^ dxJ^[hn{r,X) + hi2{r,X)]4TTr^dr (14) 

The resulting "quantum" temperatures Tq could be fitted 
to the form: 

Tq/Ep ^1.0/ {a + b^rs + cr,) (15) 

The QMC results for Ec obtained from different QM 
methods differ, e.g., by at rs = l. We used the 

most-recent Ortiz-Ballone i?cr4ata for the paramagnetic 
UEG from VMC and DMcN The small difference in 
Ec in these two simulations lead to slightly different fits. 
The fit coefHcients are, for DMC, a=1.594, b= -0.3160 
and c=0.0240, while for VMC a=1.3251, b= -0.1779 and 
c=0.0. Eight values of Vg, viz., = 1 — 6, 8, and 10 were 
used in the fit to Tq. At Vg —1 and 10, Tq/Ep goes from 
0.768 to 1.198. As 0, g{r) tends to g°{r). The UEG 
as Ts — ^- goes to a high-density fluid interacting via the 
Pauli potential. It is beyond the HNC-approximation. 

For any given rg, knowing the Tq from the 
spin-unpolarizfid case, we can obtain gij{r) and 
ExcirsjC^T)]^ at arbitrary, spin-polarization ( by solv- 
ing the coupled HNC equations. Unlike in many of the 
standard theories of electron fluids, the PDFs obtained 
from the HNC-procedure are guaranteed to be positive 
at all rg . In Fig. 2 we show typical results for gij (r) and 
compare them with those of QMC-simulations. Our re- 
sults are in excellent agreement with the DMC results. 
The reported gi2ir) do not go to zero for r=0. The de- 
pletion hole for gi2{r = 0) from DMC are in agreement 
with ours and are deeper than those from VMC.Efl The 
difference in VMC and DMC, even at rg —1, Fig. 2(a), is 
a warning that the even when Ec are in close agreement, 
other properties may have significant errors. In Fig. 2(c) 
we show the paramagnetic distribution functions at rg 
=1 and 10 together with those of the DMC simulations. 

At present there is no reliable finite-temperature mi- 
croscopic theory to compare with the g{r,T) obtained 
by our HjNTC approach. The theory of Tanaka and 
Ichimaru,Q does better than STLS and is comprehensive. 
However, like in STLS, the g{r) becomes negative for 
some values of r even at = 5. Figure 2 (d) give the 
g(r) at rg 5, T ^2Tf obtained from HNC. Unlike in SLTS 



or Tanaka and Ichimaru, the gij(r) is positive definite, as 
expected. Finite-temperature systems will be discussed 
more fully in a future publication as theji-.are relevant to 
doped-semiconductors and hot plasmas .0 

The Tq determined from the unpolarizcd Ec is used to 
calculate E^drg, C,T) at arbitrary (. The QMC results 
for £'xc(?'s,C) T=0 agree with ours, as expected from 
the agreement of our gij{r) with those from MC. For 
example, at rg = 10, the spin-polarized Ec from Ceperly- 
Alder is -0.0209 Ry, our HNC procedure gives -0.0211 
Ry, Lantto-FHNC gives -0.m|6, while KaUio and Piilo 
report a value of -0.0171 RyM 

In the theory of the UEG the response x{k,uj) of the 
interacting system is written in terms of the noninter- 
acting response function x'^{k,uj) and a LFC denoted by 
G{k,Lo). Considerable effort has been spent in obtaining 
even the static approximation G{k) for the LFC|-|Thus 
the main thrust of STLS,lI Ichimaru and Utsumijj Gel- 
dart and Taylor,B and others has been to provide the 
G{k) as a function of rg. In the static case, for a one- 
component fluid we have: 

- VuG{k) = 14-1- l/x(fc) - l/x"(fc) (16) 

where Vk = 47r/fc^. For a classical fluid we also have 

S.,{k) = -{l|P)x^,{k)/{n,n,Y'^ (17) 

Hence the LFCs can be expressed in terms of Sij [k) and 
compactly in terms of the short-ranged k-space direct 
correlation functions Cij{k). For the paramagnetic case: 

- VkG{k) = {5n(fc) + 5i2(fc) - £;i(fc)}/2/3. (18) 

We display the T=0 LFC for rg — 5 and compare it with 
the results of QMC and other theories in Fig. 3. In HNC, 
and in standard SLTS, the fc = limit is overestimated 
due to the IC. The Vashista-Singwi version of SLTS, as 
well as UI have parametrizations which externally impose 
the k=0 and k ^ 00 sum rules. The latter sum rule sates 
that G{k) ^l-g(O) for large k. The G(k) from QMC does 
not seem to follow this and is close to the second-order 
perturbation calculaton of Geldart and Taylor .0 We have 
ascertained numerically that our HNC generated G(k) 
correctly recovers the k — > cx3 sum rule very well. 

In conclusion, we have presented a simple classi- 
cal mapping of a quantum Fermi liquid of arbitrary 
spin-polarization and temperature T and shown that it 
quantitatively recovers the QMC pair-distribution func- 
tions. The only parameter is a temperature mapping 
of the correlation energy. It seems to be sufficient to 
map the static properties of the UEG. Using the method, 
we have examined spin-dependent correlation energies, 
pair-correlation functions at zero and finite T, as well as 
the local-field correction to the response functions. The 
method clearly has potential applications to Bose fluids, 
2-D electron and other systems. 
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FIG. 1. The non-interacting PDF, sr?i(r), and the Pauh 
potential, Eq. 8, at T/Ef~0 and 2 (dashed lines). They are 
universal functions of rkp, where kp is the Fermi wavevector. 



FIG. 2. The interacting PDFs gii{r) and gi2{r) at rs = l . 
Solid hnes-HNC, boxes:-DMC (a), HNC, VMC (dashes), and 
DMC. Panel (b), = 5, DMC and HNC (c), (d) are for 
T/Ef=2. In (c) The paramagnetic g{r) at rs = l and 10 at 
T=0 compared with the DMC simulations. 

FIG. 3. The local-field correction G{k), Eq. to the stati« 
rsponse at T=Q and ''s = ff Results from the HNC, QMcjia 
Geldart and Taylor (GT),|| Vashista and Singwi(VS),Q, and 
Utsumi and Ichimaru (UI)p are shown. 
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